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NUTATION  damper  -  A  SIMPLE  TWO -BODY  GYROSCOPIC  SYSTEM 

ABSTRACT 

Two  rigid  bodies  are  pinned  together  on  an  axis  which  is  a  prin¬ 
cipal  axis  of  inertia  for  each  body;  they  are  free  to  spin  about  this 
axis,  except  for  some  friction  between  them;  the  system  is  free  in 
space;  one  of  these  bodies  is  inert tally  asymmetric  about  this  mutual 
axis;  and  -  for  simplicity  -  the  other  body  is  symmetric.  If  this 
assembly  Is  inertially  more  nearly  like  a  disc  than  a  spindle,  it 
settles  in  such  a  position  that  the  mutual  axis  is  aligned  with  the 
stationary  vector  of  the  angirlar  momentum  of  this  system.  Thus  it 
constitutes  a  means  of  Internal  and  passive  damping  of  the  random 
nutations . 

The  behavior  of  this  system  is  inspected  by  considering  the  numer¬ 
able  solution  of  an  example  (in  which  one  of  the  two  bodies  is  a  thin 
disc,  and  the  other  is  a  thin  rod  which  is  always  aligned  with  some 
diameter  of  the  disc ) .  The  nature  of  the  solutions  is  Indicated 
sufficiently  clearly  by  isometric  views  of  the  trajectories  of  the 
vector  of  angular  momentum  of  the  system  in  a  coordinate  system  fixed 
with  respect  to  the  asymmetric  body  (these  trajectories  are  termed, 
somewhat  arbitrarily,  "polhpdes" ) .  The  two  extremes  of  the  friction 
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(the  "clamped"  and  the  frictionless  assembly)  are  outlined,  and  it 
is  shown  that  in  a  vaxiation  of  the  friction  between  these  two  ex¬ 
tremes  there  occurs  a  rather  peculiar  change  in  the  behavior,  of  the 
system;  a  form  of  nutation  corresponding  to  the  "unstable"  spin  of 
a  rigid  body  becomes  "metastable".  The  emphasis  of  the  presentation* 
is  on  the  legitimacy  and  usefulness  of  a  rather  primitive  "empirical" 
mathematical  approach. 
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1,  INTRODUCTION 


This  paper  describes  a  nonlinear  system  of  four  ordinary  first- 
order  differential  equations  which  represent  the  behavior  of  a  rather  ■ 
basic  assembly  of  two  rigid  bodies -a  device  which  seems  to  be  of 
particular  interest  because  by  virtue  of  its  simplicity  it  may  be 
ranked  next  to  the  classical  problem  of  the  nutations  of  a  rigid  body. 

This  device  has  been  suggested  (by  Dr.  C.  J.  Cohen  of  Naval 
Weapons  Laboratory,  Dahlgren,  Virginia,  and  others),  as  a  "nutation 
damper"  for  observation  satellites.  It  constitutes  a  passive  system 
(one  containing  friction,  but  involving  no  power  sources,  no  loss  of 
mass  of  the  satellite,  and  no  external  forces  -  such  as  magnetic) 
which  converts  the  random  nutations  of  the  satellite  into  a  smooth  spin 
about  the  vector  of  the  angular  momentum  of  the  system.  As  will  be 
discussed,  such  a  sjrstem  can  stabilize  an  assembly  which  -  inertially  - 
is  roughly  of  an  oblate  shape.  An  ancient  (though  apparently  never 
specifically  formulated)  desideratum  of  ballistics  is  an  "internal 
stabilization"  of  a  spindle.  However,  it  appears  -  so  far  -  that  a 
prolate -shaped  assembly  can  only  be  "destabilized"  by  the  system  here 
considered. 


_ f _ _ _ 

*  In  their  review  of  nonlinear  mechanics  (l),  ,  Leimanis  and  Mlnorsky 
define  as  a  basic  class  of  nonlinear  problems  that  of  the  motion  of 
a  rigid  body  fixed  at  a  point.  Basically,  these  are  third-order 
problems  (e.g.,  a  gyroscope  in  gimballs  is  a  three-body  assembly, 
but  has  only  three  degrees  of  freedom) .  To  that  class  one  may  add 
an  important  sub-class,  the  theory  of  projectile  stability  in 
ballistics;  for  a  projectile  is  usually  considered  as  a  rigid  body, 
and  -  though  the  problem  is  often  linearized  -  many  complications 
arise  from  the  mathematical  nonlinearity  of  the  forces  acting  upon 
the  projectile  (2).  The  problem  is  usually  construed  as  one  of  the 
fourth-order,  but  of  a  particular,  "eplcyclic"  subclass  of  that 
order,  which  reduces  to  the  second  order  in  complex  variables  • 

(5).  On  the  other  hand,,  many  other  radical  complications  of  the 
general  ballistic  problem  arise  in  those  cases  where  the  projectile 
may  no  longer  be  considered  a  rigid  body  (as  just  one  example  of 
this  class  of  problems  we  may  mention  a  projectile  containing  liquid), 
and  even  in  the  case  where  the  projectile  may  no  longer  be  considered 
a  body  having  inertial  axial  symmetry.  For  these  reasons,  that 
large  extension  of  Leimanis -Minorsky' s  classification,  the  motion  of 
systems  of  bodies,  will  eventually  be  of  interest;  and  in  such  a 
class,  the  problem  here  considered  seems  indeed  the  simplest  one. 
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Our  approach  to  the  problem  described  here  has  beea  '’empirical"  , 
viz.,  rather  primitive:  it  has  been  simply  an  inspection  of  the 
numerical  solutions  of  this  system  with  "sample"  inputs.  This  approach 
is  handicapped  by  imperfections  of  the  available  analog  computing 
machinery;  yet,  a  general  understeinding  of  the  behavior  of  such  a 
system  does  emerge  from  such  inspection.  In  particular,  there  have 
been  found  some  unexpected  peculiarities  in  the  behavior  of  this 
system;  they  seem  to  leave  room  for  hope  that  an  eventual  "internal- 
passive"  stabilization  of  a  spindle  might  not  be  impossible. 
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2.  STATEMENT  OF  THE  PROBLEM 


Two  rigid,  bodies  are  pinned  together  on  an  axis  which  is  a 
principal  axis  of  inertia  for  each  body;  they  are  free  to  spin  about 
this  axis^  except  for  some  friction,  between  them;  the  system  is  free 
in  space;  one  of  these  bodies  is  inertially  asymmetric  about  this 
mutual  axis;  and  -  for  simplicity  -  the  other  body  is  symmetric. 

The  standard  Lagrangian  approach  requires  an  introduction  of 
variables  of  analysis  that  completely  define,  at  each  instant,  not 
only  the  angular  velocities  of  these  two  bodies,  but  also  their 
position  (viz.,  such  coordinates  as  the  Eulerian  angles  0,  0,  \lf). 

A  considerable  simplification,  however,  becomes  possible  in  our  force- 
free  case,  through  approaching  this  problem  via  a  generalization  of 
Euler’s  dynamical  equations.  The  analysis  then  yields  only  a  part 
of  the  problem:  the  polhodes  (the  paths  of  the  vectors  of  angular 
velocity  in  the  system-fixed  coordinates);  but  this  is  the  principal 
part  of  the  problem.  Given  the  polhodes  as  a  function  of  time,  the 
position  of  these  two  bodies  at  each  instant  can  be  -  if  and  when 
needed  -  computed  through  the  Euler-Poinsot  interpretation  of  the 
motion  as  a  rolling  of  the  polhode  cones  on  the  herpolhode  cones 
(the  paths  of  the  vectors  of  angular  velocity  in  a  Newtonian  coordinate 
system);  i.e.,  by  a  subsequent  quadrature.  This  is  particularly 
convenient  in  our  force -free  case,  where  the  vector  of  the  angular 
momentum  of  the  system  Is  constant'  in  the  Newtonian  coordinate  system. 
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5.  DERIVATION  OF  THE  EQUATIONS 


Let  the  right -hajided  triads  of  the  principal  axes  of  the  two 
hodies  be  1,  2,  3  and  4,  3,  6,  with  3  and  6  being  the  common  axis 
(Figure  l).  Let  the  moments  of  inertia  be  I^  ^  1^,  I^  for  one  body; 
and  (=  l),  Ig  for  the  other.  The  components  of  the  angular 

velocities  of  the  first  body  ajre  (Hy  and  of  the  second  body, 

^4->  ^5'  ^6*  Because  of  the  constraint  on  axes  3  6  the  same 

vector  in  the  1-2 -4-5  plane  is  described  by  uii  as  by  o),  ^  ;  thus 
the  angular  velocity  of  the  second  body  can  also  be  described  by 
2  6  *  Therefore,  the  behavior  of  this  system  can  be  described  by 
a  four-dimensional  vector,  viz.,  by  the  angular  velocities  oi^,  oig, 
ojg.  Let  the  friction  torque  acting  along  the  axis  3  on  the  first  body 
be  F;  this  is  some  odd  function  of  the  relative  angular  velocity  of  the 
two  bodies 


p  S  mg  -  OJj 


For  simplicity  we  shall  here  assume  the  "viscous’’  friction  F  =  f  p  ,  f 

a  positive  constant.  I.et  also  I-,  +  I  =  K, ,  I„  +  I  3  K  ,  I,  +  I^  =  K^. 

-L  X  d  d  ^  O  ^ 

The  components  of  the  angular  momentum  L  of  the  system  on  the 

moving  axes  1,  2,  3  ai"©  related  to  the  angular  velocities  by  the 

following  equations  (the  first  two  of  which  have  the  indicated  simple 

form  because  of  the  simplifying  assumption  I^  = 

L^  =  (Ij_  +  I)  01^  = 

=  (^2  +  ^^2  ""  ^2^2 


*  However,  this  vector  suffices  to  define  the  angular  momentum  of  the 
second  body  only  in  the  cases  of  the  axial  symmetry  of  this  body,  when 
any  pair  (l,2)  of  transverse  axes  are  the  principal  axes. 
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Then  the  equations  of  motion  are 


^2  '  ^fl  h^3  =  ^ 
S  '  Va  ^  ^2“!  -  ° 
V6  *  *■ 


The  first  three  equations  express  the  conservation  of  the  vector 
of  the  angular  momentum  L  in  space i  they  are  simply  a  statement  that 

the  derivative  of  this  vector,  which  in  our  nutating  system  1,  2,  5  is 

* 

L  -  L  X  is  zero.  In  this  nutating  system  this  vector  is  conserved 
only  in  magnitude ,  and  its  direction  continually  changes .  The  fourth 
equation  is  simply  Euler's  equation  for  the  axis  6,  and  its  simplicity 
is  another  result  of  the  assumption  =  I^.  These  equations  can  be 
re -written  in  such  forms  as 


(a) 


Vl  =  "^2 


^2  =  "^1 


I^g  -  -  F 


[(K^  -  Ij)  »3  -  Vg] 

[dj  -  K,)  ] 

(K^  -  K^)  +F 


F  =  (a)g  -  oi^)  f 


(b) 


“l  "  “2 


“^2  =^^1 


0),  = 


p  -  - 


^2 

-  Kj 

“5 

M 

0'> 

^1 

^1 

-^1 

“3 

^2 

>=2 

'^2 

F 

-f 

b 

^1 

-^2 

12  I. 


jj'  =  fp 

11  TBOFNICAL  LIBSASI 

F  S  ARMY  ORDNANCE 
ABERDEEN  PEOviNG  GROOND.  MB: 
ORDBG-TL 


\  ~  \ 

['S'" 

S  ■  ^3 

* 

^2  =  -^1 

['S'" 

T.._  -  T  L.l 

~3  “5 

(c) 

-  ^l'^^ 

Lg  =  -  F  =  -  fp 

J  ■=  (I,'^  +  Ig'^)  -  Ij"'-  Lj 

Each  one  of  these  forms  has  some  advantages.  Thus,  form  (a)  might 
appear  the  most  natirral;  form  (b)  is  particularly  convenient  when  one 
needs  to  consider  the  mechanism  of  the  friction  in  detail)  and  form  (c) 
is  particularly  convenient  in  visualizing  the  problem,  for  three  reasons: 
firstly,  it  allows  a  visualization  of  L^,  L^,  as  coordinates  of  a 
point  on  a  sphere*;  secondly,  the  third  eq^uatlon  is  free  of  p;  and  thirdly, 
the  “angles  of  yaw"  of  the  axes  1,  2,  3  nutating  about  the  stationary  L 
are,  of  course,  simply  the  angles  from  the  axes  1,  2,  3  to  the  vector  L 
in  the  nutating  system  (whose  directional  cosines  are  L^/L, 

The  term  “polhode"  applies,  strictly,  to  the  paths  of  the  tips  of  the 
three-dimensional  vectors  ,  and  to,  „  and  perhaps  ought  not  be 

applied  to  the  four -dimensional  vectors  u>  „  ,  ^  and  L,  „  ,  Yet  the 
paths  of  the  latter  vectors  can  readily  be  converted  to  the  true  polhodes, 
and  do  specify  the  behavior  of  the  system  sufficiently  clearly  for  our 
purposes. 


*  This  is  particularly  convenient  in  the  work  with  an  analog  machine, 
since  it  facilitates  keeping  track  of  the  imperfections  of  the  machine. 
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4.  FORMUIATIOET  OF  THE  NUMERICAL  EXAMPLE 

Our  system  is  seen  to  be  defined  by  no  less  than  ten  parameters. 
With  our  empirical  approach  (and  for  our  descriptive  purposes)  this  is 
entirely  too  many,  and  we  need  to  simplify  the  problem  drastically,  yet 
attempting  to  lose  as  little  of  the  generality  as  possible. 

Five  of  these  parameters  characterize  the  system  "inertially" :  in 
the  forms  (a)  and  (b)  of  our  equations  they  are  ratios  such  as  Ig/K^, 

(K^  -  etc.j  in  the  form  (c),  expressions  such  as  (l/l^  +  l/Kg), 

etc.  These  five  parameters  are  restricted  somewhat  by  the  fact  that  the 
I*s  and  the  K's  are  moments  of  inertia;  these  restrictions  amount  to  the 
fact  that  the  I's  must  be  positive  and  such  that  they  can  form  a 
triangle,  viz., 

>  I  I^  -  I2I,  etc. 

It  can  be  shown  that  the  number  of  the  independent  inertial 
parameters  can  be  reduced  from  5  to  5*  However,  for  our  present  pur¬ 
poses  such  complete  systematization  is  not  imperative,  and  we  shall 

* 

not  attempt  it  here  .  Rather,  we  shall  limit  ourselves  to  a  single 
numerical  example  where 

Ig/l  =2  (i.e.,  the  second  body  is  a  "thin  disc"),  and 

Ij^  =  0,  ^2!^  ~  ^3/^  ~  (t-e.,  the  first  body  is  a  small 

"thin  rod"  along  the  axis  l). 

The  sixth  parameter  is  the  friction  parameter,  f/lg^  ^ 

in  our  three  forms  of  the  equations .  It  completes  the  specification 
of  the  system  Itself. 

The  remaining  four  parameters  are  the  initial  conditions;  that 
is,,  (in  our  three  forms  of  equations)  the  initial  values  of  0)^^  2  56’ 

“1,2,5  5,6- 


*  Some  remarks  on  departure  from  these  assumptions  are  given  in.  section 
12.  For  the  important  case  f  =  0  a  more  thorough  systematization  of 
the  inertial  parameters  of  this  system  has  been  made  by  Masaitls  (6). 
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To  simplify  our  arithmetic  we  may  select  more  convenient 
units  for  moments  of  inertia,  angular  momentum,  time  and  friction 
coefficient.  Although  distinctly  non-linear,  our  system  does  possess 
a  modicima  of  homeogeneity.  E.g.,  in  the  form  (c)  the  expressions  for 

the  derivatives  with  respect  to  time  are  linear  and  homogeneous  in  the 
reciprocals  of  the  moments  of  inertia;  hence  we  may  use  an  arbitrary 

unit  of  the  moments  of  inertia  if  we  adjust  the  unit  of  time  .accordingly. 
Furthermore,  the  first  three  of  these  expressions  are  of  the  second 
order  in  the  L's,  and  the  fourth  is  linear  in  both  L  and  f.'  As  far  as 
the  second-order  terms  are  concerned,  a  further  change  in  the  unit  of 
time  is  equivalent  to  multiplying  all  initial  conditions  by  a  constant; 
and  as  far  as  the  terms  containing  f  are  concerned,  such  change  is 
equivalent  to  a  change  of  the  friction  coefficient**  Thus  we  are  free 
to  take  an  arbitrary  unit  of  time,  provided  we  change  the  unit  of  f 
accordingly.  We  shall  select  the  units  of  time,  angular  momenta,  moments 
of  Inertia  and  f  In  such  a  way  that 

1=1  and  =  1 

We  shall  also  restrict  ourselves  to  the  natural  case  p(o)  =  0, 

which  corresponds  to  the  situation  where  the  two  bodies  are  originally 

clamped  together,  and  the  operation  of  our  mecheinlsm  starts  at  the 

** 

instant  they  are  undamped 


*  That  is,  groups  of  terms  of  the  type  dL,/dt  =  LpL,/l  can  be  inter¬ 
preted  as  dLj^/d(t/a)  =  L^L^/Cl/a),  and  further,  as  d(bL  )^(t/b)  = 

where  a  ana  b  ore  arbitrary;  while  groups  of  the  type 

dLg/dt  =  fL/l  can  be  Interpreted  as  d(bLg)/d(t/b)  =  (bf)(bL)/l,  etc. 

W^th  such  transformations  the  solutions  based  on  the  assumption 
L  =1=1  can  always  be  interpreted  in  the  conventional  units  for 
L  and  I  (the  unit  of  time  being  ab;  of  I,  a  =  I;  and  of  L,  l/b). 


**  E.g.,  the  two  bodies  forming  a  satellite  are  clamped  by  the  acceler¬ 
ation  of  the  propelling  rocket,  and  the  operation  of  the  mechanism 
starts  when  the  propulsion  ceases. 
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5.  COKSTAMTS  OF  MOTION 

f 

Two  physical  quantities  are  of  particular  interest  in  our  study. 

One  of  these  is  the  magnitude  of  the  angular  momentum  of  the 
system,-  L,  given  by 


Since  there  are  no  external  forces  acting  on  the  system,  this  is 

•H* 

a  constant  of  motion  .  This  is  one  of  the  reasons  for  our  preference 
of  the  form  (c)  of  our  equations. 

The  other  quantity  is  the  kinetic  energy  T  of  the  system,  given 
by 

2T  - 

=■  hX  ^  SX  +  (S  ■ 

where  =  I 

Generally,  this  is  not  a  constant  of  motlonj  for  it  is  readily 
seen  (particularly  from  physical  considerations)  that  T  is  in  general 
dissipated  in  friction,  at  the  rate  which  is  the  work  consumed  in 
friction.  In  fact,  by  multiplying  the  first  four  equations  of  (a) 
by  u>j^,  uOg,  tD^,  o)^,  and  adding,  we  obtain 

T=_pF=-fp^  g  0 

*  This  can  be  readily  checked  by  multiplying  the  first  three  equations 
of  (c)  by  L^,  Lg,  L^,  adding  and  integrating. 
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In  fact,  this  is  the  principle  of  operation  of  this  mechanism. 

If  there  is  any  energy  flow  because  of  the  friction,  the  system  settles 

in  a  state  in  which  p  =  0  and  T  is  a  minimum.  Generally,  one  should 

expect  this  to  occur  about  that  axis  about  which  the  moment  of  inertia 
* 

is  a  maximum  .  It  is  our  object  to  inspect  whether  such  damping  of 
the  nutation  actually  occurs,  and  if  so,  to  find  how  rapid  it  is,  and 
how  it  can  be  expedited. 

In  two  important  extreme  cases,  however,  T  is  Indeed  a  constant 
of  motion.  The  obvious  one  of  these  is  the  case  of  infinite  friction, 
when  the  system  degenerates  into  a  single  rigid  body.  The  other  is 
the  more  interesting  (and  perhaps  novel)  case  of  zero  friction  . 

Let  us  now  inspect  the  effect  of  the  parameter  f. 


*  The  underlying  principle  of  mechanics  is  invoked  often  enough  in 
specific  applications,  but  -  to  our  knowledge  -  has  never  been 
stated  in  the  full  force  and  generality  it  deserves;  it  seems  akin 
to  Schwartz's  inequality  and  to  the  second  law  of  thermodynamics. 
Cf.  also  (5). 

2 

**  In  this  case,  of  course,  any  combination  of  T,  L  and  Lg  can  be 
a  constant  of  motion.  However,  only  L  would  survive  as  such  when 
the  friction  is  introduced. 
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6.  VERY  LARGE  FRICTION 


With  f  oo  it  is  readily  seen  that  {if  the  torque  fp  is  to 
retain  its  physical  significance,  i.e.,  remain  finite)  p  -»0,  and 
the  whole  assembly  becomes  a  rigid  body  with  the  moments  of  inertia 
Kg,  andK^i  equations  (c)  become,  accordingly 

=  LgL^  (K^"^  -  Kg"^) 

Lg  =  {K^-^  -  K^“^) 

==  L^Lg  (Kg~^  -  K^"^) 

* 

The  solutions  of  these  equations  are  well  known  ,  but  will  be 
reviewed  briefly,  in  order  to  "blend”  them  with  those  of  the  more 
complicated  cases  of  0  <  f  K^oo.  The  kinetic  energy  of  the  system 
is 

T  =  L^^/2K^  +  Lg^/2K2  +  L^/2K^ 

*  The  formal  solutions  are  in  terms  of  elliptic  functions  -  which,  in 
fact,  can  best  be  defined  (4)  as  suitably-scaled  solutions  of  these 
equations.  These  solutions  can  be  summarized  as  follows;  The  sep- 
aratrlx  between  the  two  groups  of  loops  of  Figure  2  lies  in  a  plane 

LyL^  =  ^y^K^KyCKg  -  K^)/(K^  -  =  S,  say.  For  nutations  about 

axis  1,  and  initial  conditions  in  the  l-J  plane,  the  solutions  axe: 

=  L^cn(k,u)j  Lg  =  Lg^^sn(k,u)j  L^  =  Lj^Qdn(k,u).  The  modulus 

k  =  (L,  /L^)  S  is  0  at  axis  1  and  1  on  the  separatrlx.  The  inde¬ 
pendent  variable  is  u  *=  t/U,  the  unit  of  time  being  U  =  times 

/(K  -  K„)(K  -  K  )/l^.  For  nutations  about  axis  5  one  may 

V  ^  ^  d  j  J-  yj 

either  consider  k  >  1,  or  (as  is  more  customary)  take  l/k  as  modulus, 
with  L^  =  (l/k,u),  L  =  L^dn(l/k,u),  etc.  It  should  be  em¬ 

phasized  that  the  modulus  k  of  the  elliptic  functions  characterizes  not 
only  the  physical  system,  but  also  the  initial  conditions.  This  is 
somewhat  as  though  a  change  in  the  initial  conditions  changed  the 
"system")  which,  after  all,  is  the  basic  characteristic  of  non-linear 
systems . 
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and  the  tip  of  the  vector  L  (in  our  nutating  system  1,  2,  5)  therefore 
lies  on  the  surface  of  an  ellipsoid  with  seraiaxes  ^2K^  T,  , 

v/2K  T.  The  paths  of  L  are  the  intersections  of  this  ellipsoid  with 
the  sphere  L  =  const  (=  l),  and  the  motion  therefore  is  periodic. 
Figure  2  shows  these  paths  for  our  case  of  K^/'K^/K^  =  l/l.  1/2.1,  which 

.-..P  ~  .Q  4  ^  ^  />  U  .-vV.  1  4.  ^  It 


xb  cui  tsjLtmipxt;  Ui  d.  LLiaU ^bxitijjeU.^  Ui‘  UuXa.l.ti  ^  fciS36iuuly>  iu  tphiS  C8.3G 

* 

the  nutations  about  axes  1  and  5  are  stable,  while  the  spin  about  2 
is  unstable.  The  T-elllpsoid  always  has  the  same  proportions,  but  - 
it  is  very  important  to  note  -  the  magnitude  of  T  must  be  between  two 
extremes,  for  these  intersections  to  occur.  In  our  case  T  is  maximum 
for  a  pure  spin  about  1,  when  the  ellipsoid  la  wholly  outside  the 
sphere,  touching  it  only  at  its  sma.llest  semiaxis  ^2K^i  and  it  is 
minimum  for  a  pure  spin  about  5j  when  the  ellipsoid  Is  wholly  inside 


the  sphere,  touching  it  only  at  its  largest  semiaxis 

Incidentally,  from  the  condition  p  =  0  in  the  equations  (b)  it 
is  readily  seen  that  the  torque  F  necessary  to  hold  these  two  bodies 
"frozen"  is 


F  =  -  (K.  -  K-Vfl  +  I„/l.'). 

±  -d  '  ±  d’  ’  '  y  b' ' 


so  that  the  inequality  (or  I^  /  I^j,  or  (as  J.  Sternberg 
put  it)  an  "incompatibility"  of  these  two  bodies,  la  a  necessary 
condition  for  invoking  the  mechanism  of  this  friction. 


*  In  a  Newtonian  coordinate  system,  of  course,  it  is  the  axis  1 
which  nutates  about  the  stationary  Lj  etc. 
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7.  FAIRLY  LARGE  FRICTION 

With  f  finite  hut  sufficiently  large  we  should  expect  that  the  , 
motion  of  the  two-body  system  would  be  somewhat  like  that  of  a  rigid 
body,  and  the  trajectories  of  L  (or  w)  would  be  similar  to  those  of 
Figure  2j  except  that  these  trajectories  will  keep  "slipping"  towards 
the  loops  of  the  lower  kinetic  energy,  viz . ,  generally  from  axis  1  to 
axis  3*  The  assembly  thus  becomes  a  "nutation  damper",  ending  always 
with  the  pure  spin  about  axis  5.  This  is  illustrated,  for  the  case 
f  =  10,  In  Figure  3  which  has  been  drawn  by  the  differential  analyzer. 

Pour  features  distinguish  Figure  3  from  Figure  2; 

First,  the  multiplicity  of  the  trajectories  of  Figure  2  is  replaced, 
in  effect,  by  a  single  long  trajectory  (or  by  any  trajectory  intermediate 
between  any  two  swings  of  this  long  trajectory).  Of  course,  a  motion 
starting  in  a  close  proximity  to  axis  1  would  be  difficult  to  obtain  on 
a  differential  analyzer,  for  such  motion  would  be  slow  in  starting 
(it  might  best  be  obtained  by  running  the  trajectory  backwards). 

Second,  the  question  whether  the  motion  settles  on  the  positive 
pr*  the  negative  axis  3  is  a  c[uestion  of  whether  the  passage  from  what 
resembles  one  set  of  loops  on  Figure  2  to  the  other  occurs,  after 
"bouncing  away"  from  one  or  the  other  "saddle  point”;  for,  the. reader 
may ’sense  the  existence  of  two  such  saddle  points,  roughly  in  the 
vicinity  of  the  positive  and  the  negative  axis  2.  Thus,  when  starting 
with  small  nutations  about  axis  1,  the  whole  assembly  can  settle  with 
axis  3  either  parallel,  or  anti -parallel,  to  the  space-fixed  vector  L. 
Figures  3  ^'nd  4  (and  one  trajectory  on  Figure  3)  show  motion  ending  at 
the  "north  pole"  of  the  sphere.  On  each  figure  the  intersections  of 
the  long  trajectory  with  the  meridian  passing  through  axis  1  thus  mark 
a  series  of  points  starting  at  which  the  motion  goes  to  the  north  pole. 
The  same  would  hold  for  sufficiently  closely  neighboring  points  on  this 
meridian;  so  this  series  of  points  becomes  a  series  of  sub-arcs  of 
this  meridian.  In  the  region  of  loops  about  the  pole  these  sub-arcs 
blend  together;  but  in  the  region  of  loops  about  axis  1  the  situation  is 
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more  complicated.  It  can  easily  be  "interpolated"  (and  with  the 
differential  analyzer  it  can  easily  be  shown  "experimentally")  that 
for  some  trajectories  which  start  between  such  sub -arcs  in  this 
region  the  motion  would  "bounce  away"  from  a  saddle  point  near  the 
negative  axis  2,  and  would  proceed  to  the  south  pole.  Therefore  in 
this  region  this  meridian  can  be  divided  into  a  series  of  alternate 
sub-arcs  (and  the  whole  region^  into  a  series  of  alternate  stripes), 
starting  in  which  the  motion  would  go  either  to  the  north  or  to  the 
south  pole.  On  the  rest  of  the  sphere  the  situation  is  obviously 
symmetrical,  the  symmetry  being  one  of  a  rotation  through  l80°  about 
axis  2.  It  is  natural  to  Inquire  what  would  happen  to  those  two 
trajectories  which  are  the  precise  boundaries  between  these  two  " stripes" j 
one  might  expect  that  they  would  fall  into  one  or  the  other  saddle  point 
on  positive  or  negative  axis  2  and  would  "stay  there".  In  the  following 
(in  section  lO)  we  shall  touch  upon  this  matter  briefly;  but  here  we 
need  only  warn  the  reader  that  there  are  complications;  the  saddle 
points  are  not  necessarily  on  axis  2,  and  are  not  necessarily  points  of 
(even  unstable)  equilibrium.  However,  it  is  obvious  that  any  trajectory 
in  a  close  vicinity  to  these  two  "critical"  trajectories  may  spend  an 
extremely  long  time  in  a  close  vicinity  of  one  or  the  other  saddle  point 
near  axis  2. 

Third,  the  "undamplng"  of  the  nutations  about  1  proceeds  faster 
than  the  damping  of  the  nutations  about  J.  This  feature  of  our  mechanism, 
however,  may  be  affected  by  a  change  in  the  law  of  friction  employed. 

FinaJ-ly,  our  three-dimensional  representation  of  a  four -dimensional 

problem  is  necessarily  crude.  Properly,  the  representation  should 

include  some  statement  on  the  fourth  variable  (Lg,  or  p),  whose 

var'iatlon  corresponds  to  any  parameter  distinguishing  the  different 

* 

trajectories  on  Figure  2  . 

Such  curves  are  easily  obtained  on  the  differential  analyzer,  but 
are  omitted  here  for  simplicity.  The  variation  of  any  such  fourth 
variable  has  the  character  of  rapid  oscillations  superposed  upon  a 
slow  change;  but  it  is,  perhaps,  only  an  indication  of  this  slow 
change  that  is  of  interest  as  a  parameter  varying  along  our  long 
trajectory.  The  rapid  oscillations,  of  coxirse,  are  of  interest  as 
indicating  the  mechanism  by  which  the  "incompatibility"  of  the  two 
bodies  invokes  the  friction. 
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8.  MODERATE  FRICTION 


One  should  expect  that  betvreen  the  two  extremes  of  J^rlction, 
f  =  00  and  f  =  0,  there  might  exist  a  particular  value  of  friction 
at  which  the  damping  of  the  nutations  is  "maximized";  and  that  this 
is  to  he  sought  for  friction  f  <  10.  Indeed^  Figures  4  and  5, 
representing  the  cases  of  f  =  2  and  f  =  .5,  do  show  an  Increase  in 
the  rate  of  damping.  Two  features j  however,  may  he  noted.  The  gain 
in  the  rate  of  this  damping  diminishes,  so  that  this  rate  is  appar¬ 
ently  approaching  some  maocimum;  and  the  "undamping"  of  the  nutations 
about  1  becomes  much  more  pronounced  and  faster  than  the  damping  of 
the  nutations  about  5* 

In  fact,  the  passage  from  f  =  00  to  f  =  0  is  accompanied  by  a 
rather  puzzling,  almost  "qualitative",  change  in  the  behavior  of  the 
system:  there  arises  the  tendency  for  L  to  stick  for  a  long  time 
in  the  vicinity  of  axis  2  .  We  shall  presently  see  that  this 
complexity  can  be  traced  to  the  fact  that  of  the  three  pairs  of 
equilibrium  points  of  the  single-body  system  (the  points  on  the  axes 
1,  2,  3  in  Figure  2),  only  one  remains  truly  an  equilibrium  point  when 
the  friction  is  introduced.  This  is  the  pair  on  axis  3  (L^  =  Lg  =  p  =  0) . 
The  other  two  cease  to  be  equilibrium  points. 

This  phenomenon  might  be  best  inspected  by  considering  next  the 
other  (and  very  interesting)  extreme.,  the  case  of  f  =  0. 


*  This  tendency  seems  to  be  quite  distinct  from  (though  it  is 

undoubtedly  related  to)  the  fact  that  the  two  "critical."  trajectories 
on  Figures  3,  4,  5  would ^arrive  at  some  "saddle  points"  and  "stay 
there" . 
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9*  ZERO  FRICTION 

With  f  =0,  Lg  (or  oig)  becomes  a  constant  parameter,  and  our 
fourth-order  system  degenerates  into  another  third-order  system: 

h-^'a  [  (V"  -  ■'a'")  S  - 

4  =  -  -  Ki’h  I,  - 

5  =  Va 

It  may  be  noted  that  with  =  0  these  are,  in  effect,  the 

equations  for  the  nutations  of  a  rather  hypothetical  rigid  body  having 

the  moments  of  inertia  and  I,.  In  our  example  I^  «  K.  <  K  , 

L  d.  ^  ^  L  d 

and  this  hypothetical  body  is  a  thin  spindle  along  axis  3,  which  is 
considerably  different  from  our  device  made  rigid.  The  solutions,  of 
course,  are  the  elliptic  functions  cn,  sn,  dn  of  time. 

On  the  other  hand,  the  terms  in  L^,  if  they  could  be  taken  by 
themselves  (e.g.,  in  the  case  of  a  very  large  Lg),  could  be  viewed 
as  forming  a  system 

h  ■  - 

K  -  * 

which  represents  the  circular  nutations  of  another  hypothetical, 

axially-symmetrlc  (and  again  considerably  different  from  our  physical 

system)  rigid  body  with  the  moments  of  inertia  "  ^2  ~ 

K,  =  I_  +  K;  the  solutions  being  the  circular  functions  cos,  sin. -  plus 
0  5 

a  constant  (for  L^).  Thus  the  functions  L^,  L^,  for  this  zero-friction 
case  might  be  said  to  have  the  following  relation  to  the  circular  and 
elliptic  functions: 

If  the  system  x  =  f (x)  yields  elliptic  functions,  and  the  system 
£  =  s(y)  yields  circular  functions,  the  functions  L^,  L^,  are  given 
by  the  system 

L  =  f (L)  +  g(L) 
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The  shapee  of  the  curves  of  L^(t),  'L^{t) ,  L^(t)  are  Indeed  strongly 
reminiscent  of  the  elliptic  functions.  The  analytical  solution  of  this 
very  interesting  case  has  heen  achieved  hy  Dr.  ^eslovas  Masaitls  of  BRL(6). 
For  our  present  descriptive  purposes  (and  for  blending  with  the  more 
general  four-dimensional  case)  the  relationship  of  these  functions  to 
the  elliptic  functions  is,  unfortunately,  complicated  (they  are  generally 
sqiaare  roots  of  the  ratios  of  fourth-degree  polynomials  in  elliptic 
functions,  although  there  are  some  simplifications). 

There  are  three  pairs  of  "equilibrium"  points  (at  which  L  =  O): 

(a)  L,  =  L„  =  0,  L_  =  +  viz.,  on  axis  3* 

lif  3~jA 

(b)  =  0,  ^  “  ^2  ~  ^6^3  ^  ~  vizi,  in  the  plane 

2  -  3,  at  Lg  =  +  =  Lg/(l  - 

(C)  Lg  =  0,  "  ^6^3'^  ]  ^ 

1  -  3,  at  =  L^/(l  -  I^K^). 

The  simplest  way  to  determine  which  ones  of  these  are  saddle  points 
(are  unstable)  and  which  are  centers  ("stable"),  is  to  consider  the 
kinetic  energy  ellipsoid,  which  in  this  case  can  be  written  as 

l/Ai  +  +  (L^  -  =  2T  -  =  M,  say,  >  0, 

* 

so  that  the  center  qf  this  ellipsoid  is  raised  by  L^,  and  its  seinlaxes 
arevK^M,  vKgM,^IW.  In  our  example  this  ellipsoid  is  roughly  a  thin 
horizontal  "disc",  with  the  major  axis  parallel  to  axis  2.  It  can  be 
visualized  from  an  inspection  of  Figures  8.  When.M  is  minimum  this 
"disc"  is  wholly  inside  the  sphere,  touching  it  (in  the  vicinity  of  its 
major-  semiaxis  ^K^)  only  at  the  points  B;  accordingly,  these  points, 

*  There  is  no  loss  of  generality  in  reckoning  positive . 


23 


in  the  vicinity  of  the  saddle  points  of  Figure  2,  become  centers. 

With  a  large  M,  the  "disc"  in  effect  merely,  slices  off  two  (or  just 
one)  horizontal  spherical  segments;  so  that  the  points  A  are  centers^ 
too.  The  points  C  of  tangency  of  the  "disc"  and  sphere  in  the  vicinity 
of  the  intermediate  semiaxis  y of  the  "disc"  -  and  in  the  vicinity 
of  what  were  the  centers  at  axis  1  on  Figure  2  -  become  saddle  points. 

As  compared  with  Figure  2.,  Figures  8  call  for  an  additional 
parameter  because  they  no  longer  utilize' the  simplification 

Figure  2. 

The  peculiarity  of  the  behavior  of  this  physical  system  is  indeed 
due  to  this  reversal;  viz., 

the  centers  on  axis  1  in  Figure  2  become  the  saddle  points  C  in 
Figxires  8. 

the  saddle  points  on  axis  2  in  Figure  2  become  the  centers  B  in 
Figures  8;  however^ 

the  centers  on  axis  3  id  Figiire  2  remain  the  centers  A  in  Figures  8 
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10.  VERY  SMALL  FRICTION 


We  have  mentioned  that  when  the  friction  is  introduced  and  we  revert 

to  a  fourth -order  system^  the  points  A  (with  the  additional  q.ualiflcation 

that  p  =  O)  remain  equilibrium  points.  For  the  points  B  and  C,  however, 

the  requirement  that  one  or  the  other  bracket  of  equations  (c)  is  zero  - 

i.e.,  that  =  L^/{l  or  1^/(1  -  I^K^^)  -  generally 

contradicts  the  requirement  (for  a  point  of  equilibrium)  that  p  =  0,  i.e., 

that  L,  =  h,(l  +  I  /l,).  Points  B  will  remain  points  of  equilibrium  only 
y  o  0  o 

in  the  case  where 


1/(1  -  I^K^)  =  1  +  lyig,  or  =  Kg 

which  is  not  the  case  in  our  example:  l/(l  -  .l/l.l)  /  1+  •l/2,  and  which 
cannot  be  the  case  with  both  the  assumption  Ij^  =  I^  and  the  requirement 
Ig  /  ly  Thus  for  Kj  /  Kg  and  f  /  0,  strictly  speaking,  equilibrium 
points  B  and  C  do  not  exist. 


Similarly,  we  may  no  longer  use  the  argument  similar  to  that  which 
we  used  in  passing  from  the  case  of  infinite  friction  to  the  case  of 
fairly  large  friction,  for  the  quantity  M  (which  determines  the  size  of 
the  T -ellipsoid)  does  not  necessarily  have  properties  analogous  to  those 
of  T  ;  viz.,  M  as  a  function  of  t  is  not  necessarily  negative -definite . 

Nevertheless,  we  may  expect  that  with  very  small  friction  the 
trajectories  will  be  somewhat  similar  to  those  of  Figures  8.  In  this 
sense  we  shall  speak  of  points  B’  as  "quasi-centers"  for  the  loops 
similar  to  those  of  Figures  8. 

The  stability  of  points  B  so  far  (with  our  approach)  could  be 
determined  only  through  numericeil  computation.  It  can  be  described 
(for  our  example  only)  as  follows: 

1.  As  far  as  the  loops  about  the  points  are  concerned,  these 
quasi -centers  are  very  stable  asymptotically,  in  the  sense  that  the 
motion  damps  to  them  very  rapidly. 
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2.  As  fai-  as  the  location  of  the  points  b'  Is  concerned,  these 
quasi-centers  are  unstable,  with  the  ajcls  2  as  the  position  of  the 
unstable  equilibrium  of  the  points  That  is,  if  the  loops  around 
settle  with  the  quasi-center  at  the  axis  2  exactly,  the  vector 
L  will  remain  at  that  point;  but  more  generally,  during  these  oscilla 

f 

tions  the  point  B  slips  away  from  the  axis  2,  first  extremely  slowly, 
and  then  faster  and  faster. 

5.  The  slippage  of  the  points  away  from  the  axis  2  seems  to 

2 

proceed  along  a  particular  curve  on  our  sphere  of  L  . 

/ 

4,  This  slippage  of  the  points  passes  into  a  fast  nutation 
about  one  of  the  centers  A. 

5-  Nutations  about  the  centers  A  appear  asymptotically  stable. 
However,  the  rate  of  damping  of  the  nutations  about  these  centers 
becomes  extremely  small  -  so  that  it  becomes  difficult  to  distinguish 
with  the  differential  analyzer,  whether  the  centers  A  are  truly 
asymptotically  stable,  or  whether  there  might  not  exist  some 
asymptotically-stable  "small  circle"  about  the  poles  A  . 

These  features  can  be  seen  on  Figure  6  (f  =  .-1  )  and  Figure  7, 

(f  =  .02). 

These  empirical  conclusions  may,  of  course,  be  modified  by  a 
change  in  the  law  of  friction  employed. 


*  The  stability  of  points  A,  and  the  fact  that  they  are  the  only 
stable  points  have  been  recently  proven  by  D.  G.  Lewis  (7).  ' 
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11.  CRITICAL  FRICTION 


Now,  of  course,  an.  Increase  in  friction  increases  the  damping. 

It  is  not  difficult  to  find,  by  triads,  a  value  of  f  (in  our  example 

this  turns  out  to  be  approximately  f  »  .125)  with  which  the  behavior 

of  the  system  can  be  said  to  represent  most  sharply  the  border-line 

(or  the  blend)  between  the  types  of  behavior  represented  by  Figures 

>  \ 

5  and  6.  The  behavior  may  perhaps  be  described  as  follows.  Motion 

starting  in  the  "polar  zone"  of  our  sketches  passes  at  once  into  a 

damping  nutation  about  axis  3*  Motion  starting  in  the  "tropical" 

and  "temperate"  zones  rapidly  "drops"  to  the  curve  mentioned  in  item 

3  above,  and  practlcadly  stops  there  for  a  while  -  and  then  proceeds, 

with  increasing  velocity,  into  a  daxoping  nutation  about  axis  3. 
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12.  REMARKS  AND  DISCUSSION 


We  feel  that  a  differential  analyzer  -  provided  it  is  adjusted 
very  carefully  (and/or  provided  with  a  number  of  auxiliary  compensating 
circuits)  -  is  an  excellent  tool  for  the  rough  "empirical"  exploration 
of  this  sort.  Its  advantage  over  the  digital  machinery  is  the  ease 
of  access,  of  visualization,  and  of  inspection  of  the  "blends"  between 
those  types  of  behavior  of  the  system  which  at  first  glance  appear  to 
be  basically  different  from  previous  experience.  Its  shortcoming 
(in  addition  to  the  general  shortcoming  of  the  "empirical"  mathematical 
approach)  is  a  need  for  careful  adjustments,  and  an  apparent  sensitivity 
of  the  machine  to  the  type  of  circuitry  chosen. 

Our  example  was  limited  to  the  body  4,  6  being  a  "thin  disc" 

(ig  =  2l).  A  few  studies  were  done  also  with  a  "fat  disc"  (l  <  Ig/l  <  2). 
Surprisingly,  the  fat  disc  performs  much  better  than  the  thin  disc:  the 
oscillations  are  slower,  and  the  damping  is  far  more  pronounced.  As  the 
disc  approaches  a  sphere  (lg->l),  the  vector  L  appears  to  start  approach¬ 
ing  the  axis  3  asymptotically  and  without  nutations.  As  the  disc  becomes 
a  spindle  (Ig  <  l),  however,  the  motion  passes  into  a  slow  and  quicily- 
damplng  nutation  about  axis  2.  Thus  this  mechanism  allows  a  "stabiliza¬ 
tion"  of  a 'spindle  only  in  this  sense:  the  spindle  can  be  put  eventually 
into  a  pure  cartwheeling  motion,  in  which  the  "inertial  member"  (the 
small  body  1,  2,  5^  which  is  long  along  axis  1)  arranges  Itself  perpen¬ 
dicular  to  the  final  axis  of  rotation  (the  axis  2),  but  in  which  the 
orientation  of  this  axis  with  respect  to  the  spindle  (the  principal 
body  4,  5;  6)  Is,  unfortunately,  aleatory. 

A  number  of  refinements  of  this  system  are  possible.  At  this  time 
we  visualize  only  the  following  refinements: 

(a)  A  more  thorough  inspection  of  the  inequality 

1/(1  -  I^K^)  /  1  +  lyig  , 

which  as  we  have  noticed,  is ''the  reason  for  the  instability  of  the  points 
B  (and  which  represents  an  approach  to  the  condition  =  K^). 
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(b)  An  exploratory  inspection  of  the  case  where  the  body  1,  2,  5 

is  the  principal  body  ("a  fat  elliptical  disc"),  with  the  body  4,  6 

the  small  (symmetrical)  "inertial  member". 

(c)  An  exploratory  inspection  of  the  possibility  of  reinforcing 
the  inertial  action  of  the  assembly  of  the  type  (b)  by  "flywheel  gearing" 

(d)  An  exploratory  inspection  of  the  effect  of  variation  of  the 
law  of  friction  employed. 

While  this  paper  is  limited  to  purely  passive  damping,  the  "empirical" 
approach  based  on  these  principles  should  eventually  facilitate  the  inspec¬ 
tion  of  more  complicated  mechanisms,  which  may  possess  springs,  controls, 
power  sources,  and/or  loss  of  mass  of  the  satellite. 

From  the  ballistic  viewpoint,  the  most  interesting  feature  of  this 
device  la  its  strong  tendency  to  stick  -  temporarily  -  with  the  spin 
about  the  diameter  of  the  disc  (axis  2).  It  is  this  feature  -  as  well 
as  the  possibility  of  constructing  an  equivalent  of  "negative  friction" 
(viz,  an  energy  source)  -  that  seems  to  leave  hope  that  an  eventual 
"internal"  stabilization  of  a  spindle  might  be  possible.  However, '  that 
is  a  matter  of  "synthesis",  l.e.,  of  basic  changes  in  the  physical 
system  considered  here. 

SERGE  J.  ZAROODM 


*  A  crude  experimental  model  of  such  a  "flywheel  gearing"  device  shows 
an  excellent  "nutation  deanping"  performance .  It  is  interesting  to 
note  that  the  equations  of  motion  of  such  device  are  basically  the 
same  as  considered  here;  only  the  expressions  for  the  constant 
coefficients  are  different. 


29 


REFERENCES 


1.  Lelmanis,  E.  "Sonse  Recent  Advances  in  the  Dyneimlcs  of  Rigid  Bodies 
and  Celestial  Mechanics",  Part  I  of  Vol,  II,  "Dynajailcs  and  Non-Linear 
Mechanics",  of  "Surveys  in  Applied  Mathematics",  John  Wiley  &  Sons, 

Inc.  (1958). 

2.  Zaroodny,  S,  J.  "An  Exact  Re -statement  of  the  Equations  of  Motion 
of  a  Rigid  Body  for  Use  With  Modern  Computing  Machinery",  BRL 
Memorandum  Report  856  (195^) 

5.  Murphy,  G,  H.  "An  Interesting  Fourth -Order  Differential  System", 

Am.  Math.  Monthly,  LXIV,  9,  November  1957 . 

4.  Zaroodny,  S.  J.  "An  Elementary  Introduction  to  Elliptic  Functions 
of  a  Real  Variable  Based  on  the  Theory  of  Nutation  of  a  Free  Rigid 
Body",  one  of  a  series  of  ditto’d  working  papers  on  this  problem, 
available  at  BDPB,  EBL,  BRL,  APG,  Md.j  also,  Am.  Math.  Monthly, 

Aug  -  Sep  1961. 

5.  Haseltine,  W.  R.  "Passive  Damping  of  Wobbling  Satellites,  Part  I. 
General  Stability  Theory  and  First  Example",  NAVORD  Report  6579  (1959) ♦ 

6.  Masaltls,  "On  the  Motion  of  Two  Linked  Rigid  Bodies",  BRL  Report 
1097,  March,  i960. 

7.  Lewis,  D,  C.  "Paper  Without  Title",  a  manuscript  from  RIAS,  Baltimore, 
Maryland^  January  i960. 


50 


LIST  OF  SYMBOI5 


I  moment  of  inertia  of  a  rigid  body 

K  sum  of  the  two  moments  of  inertia  about  the  same  axis 
L  euigular  momentum 

T  kinetic  energy 

t  time;  dot  indicates  d/dt 
a>  angular  velocity 

P  -  -  0)^ 

f  friction  torque  coefficient 
F  .  friction  torque 

M  in  the  frictionless  case,  2T  -  Lg  /ig 
Subscripts ; 

J 

Numbers  indicate  an  axis  of  the  coordinate  system 

A,B,C  indicate  the  equilibrium  points  in  the  frictionless  case. 
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FI6.  .I-SCHEMATIC  ILLUSTRATION  OF  THE  PRINCIPLES  OF  A 
NUTATION  DAMPER,  SHOWING  THE  COORDINATE  SYSTEM  IN  WHICH, 
ON  THE  SUBSEQUENT  FIGURES,  THE  MOTION  OF  THE  TIP  OF  THE 
VECTOR  L  OF  ANGULAR  MOMENTUM  IS  GIVEN. 
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THE  RATE  OF  "UNDAMPING"  ABOUT  AXIS  I 
HAS  INCREASED  AND  SO,  TO  A  LESSER  DEGREE,  HAS  THE 
RATE  OF  DAMPING  ABOUT  AXIS  3. 

3-5 


IMPERFECTION  OF  MACHINE: 
TRAJECTORY  APPROACHES  FROM 
OUTSIDE  AND  SPIRALS  IN. 


/\ 


FIG,  f  POLHODES  FOR  VERY  SMALL  FRICTION  (f=0.02) 


THE  QUASI- CENTER  ,  B\TS  NOW  CLEARLY  DEFINED 
AS  LOOPS  THIS  POINT,  THE  POINT  SLIPS  AWAY  FROM 
AXIS  2,  SLOWLY  AT  FIRST,  BUT  WITH  INCREASING  SPEED. 


FIG.8q;  POLHODES  FOR  Z£R0  FRiCTION  (f»O.L>»«Oj 

SYSTEM  IS  A  (HYPOTHETICALLRIGID  BODY.  T  IS 
CONSTANT.  POINTS  A  AND  B  ARE  CENTERS,  POINT  C 

IS  A  SADDLE. 
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FIG.a^.  POLHODES  FOR  ZERO  FRICTION  (f  =  0.  L«  =  0.75) 
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